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ABSTRACT 

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is 
a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with 
the farthest knots. Here we present ALMA Science Verification data of CO /=2-l and 7=3-2 emission which are spatially coincident 
with the blue shifted jet of HH knots, and off'set from the disk by -18.6 km s~^ The emission has a double corkscrew morphology and 
extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew 
as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet 
which is moving at much higher velocities. We show that the 7=3-2 emission is likely heavily filtered by the interferometer, but the 
7=2-1 emission suff'ers less due to the larger beam and sensitivity to larger scale structures. Excitation analysis suggests temperatures 
exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10~^ M©, 10""^ M© km s~^ and 
^q40 gj-g respectively The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star. 
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!1. Introduction 

Disks are crucial for the formation of stars because they are a 
mass reservoir for building up the star and any planets forming 
in the system. They mediate accretion onto the star, and both 
Jets and winds from the disk act as release mechanisms for an- 
gular momentum during accretion. Jets and winds are generally 
seen qui te early on in the star formation process (e.g Class 
sources. I Arce et al]|2007l) . but have also been observed in the 
final stages of the fo rmation process as well (e.g. TW Hydra, 
iPascucci et al.ll20llh . For a complete theory of star formation, 
.we need to understand the mechanisms responsible for jet and 
wind launching. We need to understand how much of the disk 
material ends up on the star, and how much is removed from 
the system to constrain how efficiently the star is gaining mass. 
There is already evidence for residual rotation in some outflows 
launched from disks (e.g. Launhardt et al. 2009), and when the 
kinematics of the jet and disk are both measured, their rotations 
appear to be consisten t with each other (e.g. DG Tau. iTesti et al.l 
|2002; Ba cciotti et al.i r2QQ2). showing the causal links between 
these phenomena. HD 163296 is a good candidate for studying 
this phenomena because of the favourable viewing angle and the 
presence of the HH 409 system. 



HD 163296 is an isolated, prototypical Herbig Ae star which 
shows an infrared excess and is surrou nded by a Keplerian disk 
(llseUa et al.ll2Q07l : lHughes et al.ll2008h. It is at a distance o f 122 



pc, and has an age of 4 Myr ("van den Ancker et al."l998'). The 
scattered light observations of Wisniewski et al. (2008) suggest 
time variability in the disk, which they attribute to inflation of 
the inner disk edge, consist ent with the infrared variability de- 
scribed in lSitko et al.l (l2QQ8l) . These authors also show evidence 
for an emission flare at 3 jim in 2002, which they suggest may 
have been due to a combination of inner disk flare-up and ma- 
terial bein g released from the di sk into the surroundings. Using 
FUV data jDeleuil et al.l (l2005h deriv ed a mass outflow ra te of 
7 X 10"^ Mq yr"^ from the star, while iHubrig et al.l (l2009l) and 
iGarcia Lop ez et al? (zOOo) find mass accretion rates of ~ 10"^ 
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Mo yr-^ [Kraus et al. (2008) presented VLT/AMBER resolved 
Bry observations of HD 163296 and suggest that the extended 
Bry emission comes f rom a stellar wind or disk wind. The NIR- 
SPEC observations of ISalvk et all (l201ll) show rovibrational CO 
emission knots at approximately ±100 km s"^ at a velocity res- 
olution of 12.5 km s"^ which is not included in their disk model, 
but may be evidence of a high velocity wind. 

HD 1 63296 is a lso associated with Herbig Haro (HH) ob- 
ject 409 (iDevine et al. 2000). HH 409 is a chain of HH knots 
traced by [Sii] and Ha seen in both a jet and count er jet which 
is oriented perpendicular to the disk plane (e.g. .Wassell et al.l 
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I2QQ6I) . ISitko et alJ (l2008h suggested that these knots are ejected 
on 5-6 yr interv als. The PS F subtracted coronographic images 
of the region by lWassell et al. (2006) show a series of HH knots 
which they labelled A through D. The jet is blue shifted, and the 
knots are labelled A and A2. The counter-jet is redshifted, and 
the knots are labelled B-D. For the brightest components (A in 
the jet, and C in the counter jet), multi-epoch observations were 
used to derive proper motions of 0/'49 yr"^ for knot A, and 0^35 
yr"^ for C, with jet velocities greater than 250 km s"^. No sub- 
mm molecular outflow counterpart to this jet has been imaged 
before howev er there m ay be evidence for a blue shifted wind in 
the spectra of iThi et alJ (12004). 

To clarify our general terminology, we refer to a jet as high 
velocity gas (> 100 km s"^) launched from the in ner 0.1 au 
of the disk, coUimated by the magnetic field (e.g. ' Ray et afl 
120071 : iReipurth & BaUvl l200ll) . Jets are generally observed 
in optical/near-infrared atomic lines, and bow shocks resulting 
from the interaction of the jet with cloud material are often ob- 
served as HH objects (^e.g. lReipurth & Ballvll200lh . We refer to 
outflows as the molecular gas entrained by jets and/or winds as 
they push through the ambient material (e.g. Arce et al. 2007). 
These flows are generally observed to have much lower veloci- 
ties (< 50 km s"^) because they comprise ambient mater ial which 
has been swept up by the jet (JKonigl & Pudritzl [20001) . Within 
this framework, we define disk winds as the low energy material 
launched from the disk at radii of a few au (e.g. Pudritz et al. 
|2007). A disk wind (< 25 km s"^) is a centrifugally driven, mag- 
netically coUim ated wind blown off' the surface of the disk ( Pu- 
dritz & Norman Tl983l) . It is distinct from an outflow, and can 
indeed entrain an outfl ow of its o wn. An overview of these defi- 
nitions is given in Mc Kee & Qstr iker (2007). 

Here we present ALMA Science Verification data of two CO 
transitions, /=2-l (Band 6, 230.538 GHz) and /=3-2 (Band 7, 
345.796 GHz) of the disk wind associated with HD 163296 and 
HH 409. The disk gas emission is presented separately (Math- 
ews et al., de Gregorio-Monsalvo et al., submitted, and Chapil- 
lon et al., in prep) and will not be discussed here. In Section O 
we summarize our observations and data reduction. In Section [3] 
we describe the morphology and energetics of the observed disk 
wind. In Section |4] we relate these findings to the disk and jet 
properties previously determined for this source, and summarize 
in Section [5] 



2. Observations and Data Reduction 

These ^^CO /=2-l and /=3-2 observations were taken 
with ALMA as part of the Science Verification program 
(2011.0.000010.SV). Both datasets were observed in June and 
July 2012 with 18-20 antennas in an extended configuration 
(-20-400 m baselines). The Band 6 data consisted of three sep- 
arate observing blocks, and the Band 7 data are a combination 
of five. The total time on source in each band was 46 and 59 
minutes, respectively. The quasars J 1924-292 and J1733-130 
were used for bandpass and phase calibration. For the Band 6 
data, Juno, Neptune and Mars were used for amplitude calibra- 
tion. For four of the Band 7 executions, Neptune was observed 
for amplitude calibration, and this amplitude calibration was ap- 
plied (bootstrapped) to the fifth execution for which Juno was 
observed, but was too close to the horizon to be properly used for 
amplitude calibration. Data reduction and imaging were all per- 
formed using the Common Astronomy Software Applications 
(CASA, McMullin et al. 2007). Both data sets were cleaned us- 
ing natural weighting, and the channel spacings of the spectral 
windows for the CO lines were 244 kHz and 122 kHz, which 
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Fig. 1. Visibility spectra averaging over the two shortest baselines in 
each of the Band 6 and Band 7 datasets showing the blue and red shifted 
wind emiss ion in this field of v iew. The rest velocity of the source is 
5.8 km s~^ (iHughes et al.ll2008h as shown with a green line in the plot. 
The Band 6 (CO 7=2-1) data have been shifted upwards by 5 to clearly 
show both spectra. 

Table 1. Properties of the 230 GHz and 345 GHz Observations. The 
rms noise values listed below have units of mJy beam~^ per channel. 
The channel spacings are 0.32 and 0.25 km s~^ for the 230 and 345 
GHz data (respectively) 



Line 



rms 



Synthesized 
Beam, P.A. 



Primary 
Beam 



Largest 
Observable Scale 



/=2-l 
/=3-2 



6.9 
4.2 



0.80' 
0.68' 



xO.68'', 75° 
xO.44'', 91° 



26" 

17'' 



11' 



were resampled to spectral resolutions of 0.32 and 0.25 km s"^ 
in Bands 6 and 7 respectively. Both datasets were self-calibrated 
in both phase and amplitude using line-free continuum channels 
in all four spectral windows. The per-channel rms noises in these 
two datasets are 6.9 and 4.2 mJy beam"^ after self-calibration. 
Data from two antennas (PMOl and DV03) were removed for 
the July observations in Band 6 because of anomalously high 
amplitudes. HCO^ /=4-3 was simultaneously observed with CO 
J=3-2 in the Band 7 observations, and thus has the same noise 
properties. These data are presented in Mathews et al. (Submit- 
ted). 

The synthesized beams for the two sets of observations are 
shown in Table [T] The shortest baselines in the array dictate the 
largest scale observable in each dataset, and the formula for de- 
termining the approximate largest angular scales (LAS) is siven 
in Section 6.2 of the ALMA Cycle 1 Technical Handboolfl The 
largest scale structures to which our observations are sensitive 
are given in Table \T\ Any structures larger than this are filtered 
out by the interferometer. 



3. Results 

Figure [T] shows the averaged visibility spectra of the two short- 
est baselines (21.2 and 21.7 m) in both datasets. The velocity 
integration limits used for finding emission in the image plane 



^ https://almascience.nrao.edu/documents-and-tools/cycle-l/alma- 
technical-handbook 
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Fig. 2. Blue shifted disk wind from HD 163296. The greyscale and single (4cr) contour show the CO /=2-l integrated intensity of the blue 
shifted wind from HD 163296. The blue and red contours show the 15, 20 and 25 times the rms noise (20 mJy beam~^) of the blue (-2 to 5.5 
kms~^) and red (6 to 12 km s~^) HCO^ 7=4-3 emission in the disk. The inset image shows the first moment map of the HCO^ 7=4-3 emission 
in the disk, with the same red and blue contours as in the main figure for reference. The two black arcs show the primary beams of the 345 GHz 
(smaller) and 230 GHz (larger) observations. The cyan contours show 4, 6, and 8 times the rms noise in the CO 7=3-2 emission (2.6 mJy beam~^). 
The two solid purple circles show the 2012 positions of knots A2 and A3 as described in the text. The dashed purple circle and arrow show the 
position of knot A2 in 2004, and its projected motion to 2012. The two dark cyan circles in the bottom left comer show the synthesized beams of 
the Band 7 (left) and Band 6 (right) datasets. The dimensions of the synthesized beams are given in Table [T] 



are shown as red and blue boxes in Figure [T] This wind femis- 
sion is at distinct velocities from the disk itself (centered iit the 
green line in Figure [TJ. The redshifted emission is strongest in 
the Band 6 data. At both the red and blue shifted velocities high- 
lighted in Figure[T] there is phase coherence, and it is offset from 
zero. This suggests that the emission is coming from a region 
offset from the observing center. 

We detect blue shifted wind emission in the image plane in 
both CO 7=2-1 and 7=3-2. Figure [2] shows the integrated inten- 
sities of the CO emission as traced by CO 7=2-1 in the gre} scale 
and in CO 7=3-2 in cyan contours. The red and blue contours in 
the center of the figure and in the inset show the 15, 20, and 25 
cr contours of the integrated blue and red shifted HCO^ j =4-3 
emission in the disk. This is shown to give context to the posi- 
tion of the wind. We have plotted (using dashed lines) the 2004 
position of knot A2 and its proper motion. Also shown (with 
solid purple circles) are the expected positions of knots A2 and 
A3 in 2012 using the proper motion of knot A given in .Was- 



sell et al. (l2QQ6l) . This correction for proper motion is described 
more ful y in Section IH The pointing center is ck= 17:56:21. 28, 
^=-21:57:22.36, and the large black circles show the full width 
at half maximum power of the Band 6 and 7 primary beams (the 
FWHM J)f the Band 6 primary beam is larger than the Band 7). 
Neither of the datasets shown in Figures |2] and [3] are primary 
beam coi'rected. The integrated and peak intensities of the CO 
7=2-1 emission are 2.58 Jy km s"^ and 0.084 Jy beam"^ km s"^ 
respectively. The disk wind traced by CO 7=2-1 is a continuous 
structure with what appears to be a forked end. From the first 
moment maps of both transitions (the CO 7=2-1 first moment 
map is sliown in Figure [3]), we find that both transitions show 
velocity gradients perpendicular to the wind and jet direction, 
consister t with that in the disk (blue-shifted to the south-east, 
red-shift(id to the north-west). 

The !>patial extent of the 7=3-2 emission is much smaller 
than the 7=2-1 emission (see Figure[2]). It is concentrated in two 
knots: on e within the full width at half maximum of the primary 
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beam, and one external to it. The brightest CO J =3-2 knot is 
well within the primary beam, and has integrated and peak in- 
tensities of 0.6 Jy km s"^ 0.06 Jy beam"^ km s"^ respectively. 
We suggest it is associated with a jet knot ejected in 2002. The 
second knot(A2), outside the Band 7 primary beam, is coinci- 
dent with where H H knot A2 should be in 2012 using the proper 
motions derived in IWassell et al.l (l2006l) . This knot is best seen 
when the Band 7 data is primary beam corrected (as shown in 
Figure 13]). Primary beam correction takes into account the atten- 
uation of the signal away from the pointing center. The noise in 
the map is no longer uniform, but increases with distance from 
center. The smaller and larger arcs in Figure 51 show where the 
Band 7 and Band 6 primary beams are attenuated to 50%. The 
primary beam corrected data extends to where the primary beam 
is attenuated to 20% of its sensitivity at the pointing center. The 
CO J=3-2 emission in the vicinity of HH knot A2 is much more 
prominent in this representation. 

The differing morphology between the two CO transitions 
may be explained by spatial filtering by the interferometer. The 
Band 7 observations cannot detect structures larger than about 
8'', while the Band 6 observations are sensitive to larger scale 
structures (~ IT')- It is a combination of the higher excitation 
conditions required for the J =3-2 transition and the spatial fil- 
tering which is causing the wind to appear clumpy and compact 
in the Band 7 observations. We masked the shortest baselines 
in the Band 6 data to match the largest angular scales recover- 
able in the Band 7 data (8'' at 230.5 GHz corresponds to 16.6 
k/1). In this filtered data, we find that the wind, as traced by CO 
/=2-l, appears clumpier, and peaks close to the J =3-2 emission. 
That the regions with the strongest emission for two species are 
slightly diff'erent suggests that the excitation of the two lines also 
plays a role in the morphology. 

In the image plane (see FigureO, we do not detect the north- 
em redshifted gas. It is likely that the redshifted wind is filtered 
out by the interferometer. This is evidenced by the redshifted 
emission being detected in the Band 6 visibilities, but not the 
Band 7 visibilities. The redshifted emission is filtered out on 
all but the shortest two baselines (at ~ 21 m). The next short- 
est baseline is 27.5 m, and the short baseline sampling appears 
to be too coarse to image the red lobe. Followup observations 
covering a wider areas, and recovering larger scale structures are 
required to confirm this hypothesis. 

4. Discussion 



Wassell et al. ('2006') characterized the knots of HH 409 associ- 
ated with HD 163296. For knot A they suggest a proper motion 
of 0.'''49 yr"^. Using that proper motion, and extrapolating to 
2012, we find that knot A is well beyond the primary beam of 
our observations. However, if we apply the proper motion of 
knot A to knot A2, we find that it should have been approx- 
imately IT' from the disk when our observ ations were taken . 
Using a position angle of 42° (taken from Wassell et al.l l2006'), 
this knot should be slightly ahead of the CO emission seen in 
the primary beam corrected image of Figure 5] as shown with a 
purple circle. Note that this clump, at 1 T', is in the region where 
the interferometer is only 20% as sensitive as it is at the pointing 
center. That a jet knot should be co-located with this CO J=3- 
2 emission feature suggests that the molecular emission is real 
despite being far from the pointing center. 

The observations of S itko et al.l (i2008j showed a flare in 3 
//m emission coming from HD 163296 in 2002. If we assume 
that this was a knot launching event, and assume this knot has 
the same proper motion as HH knot A (0.''49 yr"^), then, at the 



epoch of our observations (10 years later), the ejecta should be 
4.9'' from the disk. This is the location of the leading edge of our 
strongest CO J=3-2 knot, which we named A3 (see Figure O. 
It should be noted that an atomic HH knot from this flare would 
h ave been obscured by the coronograph of the 2004 observations 
of lWassell et al.l (l2006l) . but that subsequent obser vations should 
show the atomic knot emission (see for instance, iGUnther et al.l 
I2013L and Ellerbroek et al. in prep.). 

Based on the emission morphologies shown in Figures |2]and 
O we suggest that the ^^CO emission is primarily tracing the disk 
wind coming from HD 163296. The morphology of the emis- 
sion, and its velocity structure suggests that this is a disk wind, 
and not an entrained outflow from the jet. That there appears 
to be two arms to the wind argues against entrainment because 
it is unlikely that a single knot can produce two velocity com- 
ponents (like those highlighted by dashed lines in Figure [3]). If 
the CO were entrained material, there is no explanation for the 
highly blue shifted material at the end of the bluest arm of the 
observable wind, since there is no HH knot ahead of it which 
could have entrained that gas. In the area around the knots, the 
molecular gas in the wind is being heated by the jet; either from 
the warm shocks, or directly from the jet heating the dust. Be- 
cause the temperature is elevated, smaller amounts of gas appear 
brighter. This means that compact and bright knots of CO will 
appear in the larger scale wind. In the case of CO /=3-2, the 
larger scale wind is resolved out by the interferometer, but the 
bright and compact knots of heated gas are still detectable. The 
CO /=2-l emission is also quite bright in this region. Note that 
the enhanced CO emission is slightly closer to the disk ( i.e. 
upstream) than the expected position of the knots in 2012, rein- 
forcing the suggestion that this material was heated by a likely 
HH knot. 

The velocity structure shown in the CO /=2-l first moment 
map (FigureO suggests that there may be two foot points for the 
disk wind, which are on opposite sides of the disk, co-rotating 
with it. The velocity field is an imprint of the rotation of the disk, 
and the wind itself may be collimated by the magnetic field. The 
bluest portion of the disk wind is emitting from the south east 
part of the disk, which itself is known to be blue shifted. The red- 
dest part of the disk wind is coming from the north west portion 
of the disk, which is red shifted. Some jet s do show evidence fo r 
rotation like that seen here (e.g. DG Tau lBacciotti et al.ll2002l) . 
We suggest this disk wind has a double corkscrew, similar to the 
m odelled disk wind p henomena se en in, for instance. Figu re 3 
of lOuved et all (120031) or Figure 7 of iLaunhardt et al.l (l2009l) . but 
with two components. The dashed blue and red lines in Figure [3] 
give an illustration of the double corkscrew described here. The 
two components of the wind cross each other along the line of 
sight just further away from the disk than the A3 knot seen in CO 
J =3 2. This double corkscrew can explain why the outer edges 
of the /=2-l emission appear forked, and with very diff'erent ve- 
locities; they are at opposite en ds of the double cork screw. 

The recent observations of Gunther et al. (2013) show Lya 
emission associated with the HH knots, which indicates that the 
gas in these regions is h ot. Using RADEX, a no n-LTE ID radia- 
tive transfer calculator (" van der Tak et al.ll2007b , we determined 
the excitation temperature and column density necessary to pro- 
duce the ratio of line peak intensities of the /=2-l (1.195 K) and 
J =3-2 (1.523 K) transitions of CO integrated over the disk wind 
emission near knot A3. For consistency, we used the spatially 
filtered Band 6 data for comparison with the Band 7 data. Us- 
ing an ambient density o f 1.9x10^ cm~ ^ (see Table 8 of Wassell 
et al. l2006l) . we find that to match our line ratio requires a column 
density of 1 x 10^^ cm"^. Dividing this derived column density by 
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Fig. 3. CO /=2-l intensity weighted velocity map of the blue shifted 
wind from HD 163296. The purple circles show the 2012 positions of 
the knots, the blue and red contours show the HCO^ 7=4-3 emission, 
and the two dark cyan circles in the bottom left comer show the syn- 
thesized beams of the observations, as in Figure |2] The blue and red 
dashed lines delineate the double corkscrew described in the text. 
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Fig. 4. CO J -3-2 integrated intensity map of the blue shifted wind 
from HD 163296. The colourscale shows the primary beam corrected 
integrated intensity of the J -3-2 emission The black, blue and red cir- 
cl es and contours are the same as Figure 2 The molecular counterpart 
to IhH knot A2 is more prominent in this image because we have cor- 
rected for the fall-off in sensitivity as a function of distance from the 
pointing center. 



the assumed volume density recovers our observed emitting re- 
gion when we make the assumption of spherical symmetry. For 
higher column densities, the line ratio does not converge to the 
observed J=2-l/J=3-2 ratio of 0.8. For lower column densities, 
the line peak intensities are both below 1 K. At this column den- 
sity the line ratio approaches 0.8 for temperatures greater than 
900 K. Above 1000 K, the /=2-l emission has a negative exci- 
tation temperature. We can recover our line ratio and integrated 
intensities at lower temperature, however this results in a radius 
ten times greater than that observed when dividing the column 
density by the volume density. We therefore use an excitation 
temperature of 960 for further calculations. At these densities 
and temperatures, both lines are optically thin, with opacities 
near 0.1. These high temperatures are not unrealistic, especially 
since the region we are probing is being heated by a jet knot. 

Using the derived column density and assumed density, the 
depth of the emission region towards the brightest CO/=3-2 
knot (~ 300 au) is roughly the same as the width of the emission. 
This is consistent with the flow being cylindrical. To calculate 
the mass and energetics of the outflow, we used the integrated 
intensity (moment zero, e.g. Figure |2]) and intensity weighted 
velocity (moment one, e.g. Figure [3]) maps of the CO emission 
in both CO transitions. When creating these maps, the data were 
clipped at 3.5 and 4 cr in each channel when creating the zeroth 
and first moment maps (respectively). We assume a wind age of 
400 yr based on the average velocity in the wind (-18.6 km s"^) 
and its spatial extent in the Band 6 observations (~ 13'0- This 
will underestimate the age of the wind if it is more extended than 
the Band 6 primary beam. This is the age used to determine the 
wind forces, mechanical luminosities and mass loss rates quoted 
in Table [21 

To calculate the mass, we converted the integrated intensity 
of the line into a CO column density assuming local thermo- 
dynamical equilibrium, a temperature of 960 K and a Jansky to 
Kelvin conversion based on the synthesized beams of the Band 6 
(43.48 Jy/K) and Band 7 (33.78 Jy/K) data. This was converted 
to a total molecular column density using a CO abundance of 
10""^ with respect to molecular hydrogen, and a mean molecu- 
lar mass of 2.8 (see, for instance, appendix A of Kauff'mann et 



al."2008). The wind momentum was derived by determining the 
mean velocity of the gas (-18.6 km s"^), and multiplying that by 
the derived mass. We calculated the wind force by dividing the 
momentum by the wind age. Similarly, the mechanical energy in 
the wind was calculated using E= ^mv^, and the mechanical lu- 
minosity was determined by dividing the energy by the assumed 
age of the outflow. Using the column density derived from the 
RADEX calculation, the outflow mass drops by almost an order 
of magnitude. 

The angular momentum (L = mvr) of the A3 knot is 
3.7x10"^ Mo km^ s"^ where r is the radius of the CO /=3-2 
emission (r = 133 au). A mass accretion rate of ~ 10"^ M© 
yr"^ was previously determined for this source ( Garcia Lopez 
et al. ' 20061 ; iHubrig et al.l l2009l) . This is a factor of two lower 
than the mass loss rate we present in Table |2] (2x10"^ M© yr"^). 
In general, mass loss rates are expected to be about 10% that 
of th e accretion rate. Th e disk itself has a mass of order 10"^ 
M© (JTilling et al.ll2012h . which, when coupled with the mass 
loss rate we derive, suggests a disk dispersal timescale of 10^ 
yr. Given the age of the system (4 Myr), it is possible that the 
star is in the final stages of accretion, and is blowing out the last 
of its inf alling material via the final stages of a disk wind. We 
note that lKraus et al.l (l2008h suggest that the Bry accretion rate 
for this source may be contaminated by outflow emission, which 
further strengthens our argument, as the mass loss rate would 
then be higher than the mass accre tion rate. If the gas we re from 
a photo evaporative disk wind (i.e. lHoUenbach et aDl2000l) . then 
the characteristic velocity of the gas measured here (-18.6 km 
s"^) should be roughly equivalent to the sound speed. A sound 
speed of 18.6 km s"^ indicates a temperature of 10^ K, which is 
much too high for molecular gas to be present. Thus, we are not 
seeing a photo evaporative wind. 



5. Conclusions 

We have presented here the first evidence for a molecular disk 
wind coming from a Herbig Ae star. This disk wind is coincident 
with a chain of HH knots. We have quantified the energetics 
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Table 2. Outflow energetics for the two CO transitions. The first CO 
/=2-l column is derived from integrating over all of the CO 7=2-1 
emission, while the second is derived from the data that was filtered 
to sample the same spatial scales as the Band 7 data. This 'Filtered' 
column is more directly comparable to the 7=3-2 data. The quoted un- 
certainties come from the uncertainties in the fluxes, propagated through 
the equations in quadrature. Note that the quoted velocity is blue shifted 
with respect to the source velocity of 5.8 km s~^ 





CO 7=3-2 


CO 7=2-1 


CO 7=2-1 






Integrated 


Filtered 


M (10"*' Mo) 


5.7+0.06 


73.1+0.8 


24.6+0.3 


V(kms-i) 


-18.6+0.1 


-18.6+0.3 


-18.9+0.5 


P(10-5Mokms-i) 


10.7+0.1 


136.0+3.5 


46.4+4.6 


E (10^0 erg) 


1.98+0.04 


25.1+1.0 


8.7+0.5 


Wind kinematic age 


= 400 yr 




F« 


2.7+0.3 


34.0+4 


11.5+2 


L (10-3 L^) 


0.41+0.05 


5.2+0.7 


1.8+0.3 


Mou,(10-8Moyr-i) 


1.4+0.2 


18.3+2 


6.2+0.7 



The units of the wind force are 
10"^ Mo km s-i yr-i. 
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of the blue shifted wind emission, and show that the positions 
of the two CO J =3-2 emission knots are consistent with one 
of the optically detected HH objects (A2), and a younger knot 
(A3) which was likely emitted during a flare in 2002. Further 
followup ALMA observations, including wider-field mosaicing 
and incorporating larger scale emission will be required to fully 
characterize this phenomena. 

The mass being released from the system via the disk wind 
is twice the mass accr etion rates found via Bry emission ( Gar- 
cia Lopez et al. 2006). For s table, accreting syste ms, this ratio 
is generally closer to 0.1 (e.g. lHartigan et al.lll995h . The deriva- 
tion of the accretion rate using Bry is a tim e average (Muzerolle 
et al. 1998), and since the age used for the wind is kinematically 
derived, so is the mass loss rate in the wind. Because the mass 
loss rate is, on average, larger than the mass accretion rate, this 
system cannot be in a stable growth stage. 
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